uvbyCaHP CCD Photometry of Clusters. III. The Most 
Metal-Rich Open Cluster, NGC 6253 
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ABSTRACT 

CCD photometry on the intermediate-band uvbyCaUP system is presented for the old open 
cluster, NGC 6253. Despite a high level of field star contamination due to its location toward 
the galactic center, combination of the data from the multiple color indices with the core cluster 
sample derived from radial star counts leads to the identification of a set of highly probable, 
single cluster members. Photometric analysis of a select sample of 71 turnoff stars produces a 
reddening value of E{b - y) ^ 0.190 ± 0.002 (s.e.m.) or E{B ~V) ^ 0.260 ±0.003 (s.e.m.) from 
71 stars. The metallicity indices, Smi and Shk, both identify this cluster as the most metal- 
rich object studied on either system to date. Simple extrapolation of the available metallicity 
calibrations leads to [Fe/H] ranging from +0.7 to -1-0.9. Metal-rich isochrones with overshoot 
imply an age between 2.5 and 3.5 Gyr, with an apparent distance modulus between (to — M) — 
11.6 and 12.2, depending upon the isochrones used. The improvement in the fit using a-enhanced 
isochrones may indicate that the cluster [Fc/H] is closer to -1-0.4, but the photometric indices are 
distorted by an elemental distribution other than a scaled solar. The galactocentric position of 
the cluster, in conjunction with data for other clusters and Cepheids, is consistent with the inner 
disk reaching and maintaining a metallicity well above solar since the early history of the disk, 
unlike the solar neighborhood. 

Subject headings: color-magnitude diagrams — open clusters and associationsundividual (NGC 6253) 



1. INTRODUCTION 

The shared properties of distance, age, and 
composition for stars within a cluster supply justi- 
fication for their study within the context of stel- 
lar and galactic evolution. An example of how a 
well-defined and homogeneously- analyzed cluster 
sample can shed light on galactic evolution is pro- 
vided by Twarog et al. (1997), wherein analysis of 
over 70 open clusters demonstrates the existence 
of structure within the galactic abundance gradi- 
ent, structure that has been corroborated most re- 
cently through the use of Cepheids by Andrievsky 
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et al. (2002b). Despite the unusual combination 
of high precision and sample size for the collected 
data, critical regions of the galactic disk remain 
inadequately sampled. Moreover, only a modest 
fraction of the open clusters within observational 
reach have their fundamental parameters deter- 
mined to a degree of reliability that can satisfy 
the needs of Galactic studies tied to the global 
characteristics of the population. 

Though the number of broad-band photo- 
metric surveys of individual clusters with color- 
magnitude diagrams (CMD) of sufficient precision 
has increased significantly in the last five years 
with the expanded availability of larger-format 
CCD's, the interpretation of these data has been 
hampered by the absence of a comparable progress 
in determinations of the critically important prop- 
erties of reddening and metallicity, linchpins to age 
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and distance. Attempts to constrain the param- 
eters using CMD morphology and comparison to 
theoretically synthesized CMD's have had some 
success (Tosi ct al. 1991; Bonifazi ct al. 1990; 
Gozzoh et al. 1996; Marconi et al. 1997), though 
the answers drawn from different models do not 
always concur and the uncertainties in the pa- 
rameters remain larger than desirable for many 
projects. 

To improve this situation for individual clus- 
ters and to expand and test the sample delineat- 
ing the galactic structure outlined by Twarog et 
al. (1997), a program of CCD photometry of open 
clusters on the extended Stromgren intermediate- 
band {uvbyCa) and H/3 systems has been under- 
taken. Examples of the initial work illustrating the 
value and precision attainable with careful analy- 
sis have been given for the metal-rich open cluster, 
IC 4651 (Anthony- Twarog & Twarog 2000a, here- 
inafter referred to as Paper I) and the globular 
cluster, NGC 6397 (Anthony-Twarog & Twarog 
2000b, hereinafter referred to as Paper II). These 
clusters were selected to sample the extreme range 
of cluster parameters with objects that had been 
reasonably well-studied by other photometric and 
spectroscopic means. Our first program cluster 
with no intermediate-band photometric work or 
reliable high-resolution spectroscopic abundances 
is the distant open cluster, NGC 6253. 

NGC 6253 has been the focus of three broad- 
band photometric surveys to date, Bragaglia et 
al. (1997), Piatti et al. (1998), and Sagar et al. 
(2001) (hereinafter referred to as BR, PI, and SA, 
respectively). It is of interest because, as carefully 
demonstrated in BR, the cluster is located toward 
the galactic center (Z^^ = 336°, 6^^ = -6°), sig- 
nificantly older than the Hyades, and super-solar 
in metallicity, while exhibiting a richly-populated 
turnoff and giant branch. Despite the concurrence 
of opinion regarding the broad features of the clus- 
ter, examination of previous work reveals a wide 
range for the key parameters of interest. Previ- 
ously published work produces plausible redden- 
ing estimates ranging from E{B — V) = 0.20 to 
0.32, metallicity from solar abundance to almost 
triple the solar value, and an age from 2 to 5 Gyr; 
we will discuss the specifics of each property and 
its past derivation in Sec. 4. Since traditional 
linear abundance gradients imply a change in the 
mean [Fe/H] of less than 0.3 dex over a 4 kpc 



galactocentric range, definitive evaluation of any 
structure within 2 kpc of the sun requires [Fe/H] 
data of significantly higher accuracy than is cur- 
rently available for NGC 6253, for a large number 
of clusters. Moreover, a definitive metallicity is 
important in defining the value of the cluster as 
a template for younger, metal-rich, extragalactic 
populations and as a potential source of stars with 
planetary systems. 

The underlying basis for the use of intermediate- 
band photometry is simple. The brighter observa- 
tional limit imposed by the narrower bandpasses 
is, in part, compensated by the higher precision 
with which one may disentangle the stellar pa- 
rameters. DDO photometry (Twarog et al. 1997) 
reconinicincls itself because it focuses on the bright- 
est stars within a cluster, the red giants, and the 
calibration of the system is exquisitely tied to the 
spectroscopic abundances of the extensive sample 
of field giants near the sun (Twarog & Anthony- 
Twarog 1996). The weaknesses of the system are 
the small sample of giants often present in many 
poorly populated clusters observed within a field 
rich with cluster non-members and the possibility 
that post-main-sequence evolution may alter the 
atmospheric abundances of some stars. 

The extended Stromgren and II/3 systems ex- 
ploit the rich population of main sequence stars 
and have been well-calibrated for spectral types 
ranging from B through G (Crawford 1975, 1978, 
1979; Olsen 1988; Schuster & Nissen 1989), the 
temperature range typified by the turnoff in most 
open clusters. No open cluster has been identi- 
fied that exhibits a range in [Fe/H] or a difference 
in [Fe/H] between the turnoff stars and the red 
giants when the abundances are based primarily 
upon the Fe lines. The larger population of main 
sequence stars can lead to easier identification of a 
probable sample of cluster members while reduc- 
ing the standard error of the mean for the final 
cluster parameters. The multiple indices also en- 
hance the elimination of probable non-members 
from the final sample. For an example of how 
successful this approach can be for a typical open 
cluster, the reader is referred to Paper I. 

Section 2 contains the details of the CCD ob- 
servations and their reduction and transformation 
to the standard system. In Sec. 3 we discuss the 
CMD and begin the process of identifying the sam- 
ple of probable cluster members. Sec. 4 contains 
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the derivation of the fundamental cluster parame- 
ters of reddening and metallicity. In Sec. 5, these 
are combined with broad-band data to derive the 
distance and age through comparisons with theo- 
retical isochrones. Sec. 6 summarizes our conclu- 
sions in the context of our current understanding 
of the evolution of the galactic abundance gradi- 
ent. Preliminary results for NGC 6253 presented 
in Twarog et al. (2002) are superseded by this 
discussion. 

2. The Data 

2.1. Observations 

The new photometric data for NGC 6253 were 
obtained using the Cassegrain-focus CCD imager 
on the National Optical Astronomy Observatory's 
0.9-m telescope at Cerro Tololo Interamerican Ob- 
servatory. We used a Tektronix 2048 by 2048 de- 
tector at the //13.5 focus of the telescope, modu- 
lated by CTIO's 4" x 4" uvby filters and our own 
3" X 3" H/3 and Ca filters. The field size is 13.5' on 
a side. Frames were obtained in April 1999 with 
all filters except the H/3 filters; frames obtained in 
May 2000 included all bandpasses except Ca. In 
all, 11 to 22 frames were acquired of NGC 6253 for 
each filter with the following total exposure times: 
37, 63, 71 and 159 minutes for y,b,v and u, 83 
minutes for Ca, and 40 and 183 minutes for the 
H/3 wide and narrow filters. 

Standard IRAF routines were used to perform 
initial processing at the telescope; a fairly com- 
prehensive description of these steps and the pro- 
cedures by which we merge photometry based on 
profile-fitting algorithms may be found in Paper I. 

2.2. Reduction and Transformation 

The procedures outlined in Paper I yield photo- 
metric indices of very high internal precision. We 
outline here the steps taken to transform these in- 
strumental indices to the standard systems, and 
will separately note specific issues pertaining to 
the reduction of uvby indices, hk indices, and H/3 
indices following this overview. 

In addition to our program fields, we observe 
standard stars selected from the photometry cat- 
alogs of Olsen (1983, 1993, 1994) and Twarog & 
Anthony-Twarog (1995) and in open clusters on 
each photometric night. These are used to de- 



termine extinction solutions for each photomet- 
ric night for all indices other than H/3. Aper- 
ture photometry is obtained for every standard 
star and program cluster on all but the deepest 
exposure frames using aperture radii ranging from 
8 to 13 pixels centered on the star, along with 
sky measurements using annuli enclosing a compa- 
rable area. Following correction for atmospheric 
extinction and the exclusion of stars with the 
most serious crowding, average magnitudes in each 
bandpass are constructed for all stars. From the 
mean magnitudes, mean indices are constructed 
for stars in every field on each night. For a pro- 
gram cluster observed on a photometric night, the 
mean difference between the aperture photometry 
and the profile-fit-based photometry is then deter- 
mined for each index based on a selected sample 
of relatively uncrowded stars. 

In determining the transformation equations 
between the instrumental and standard system, 
our preferred approach is to use the photometry 
within the open clusters to establish the slopes of 
the transformation relations and, if possible, ap- 
peal to photometric catalog values of single-star 
standards to establish the zero point of the calibra- 
tion equations. With these transformation equa- 
tions, indices constructed from all the frames may 
be converted to standard values. By inference, the 
profile-fit photometry may be calibrated with the 
same relations. 

H/3 frames of NGC 6253, as well as several stan- 
dard fields in open clusters, were obtained on the 
night of May 2, 2000. The transformation equa- 
tion for conversion of aperture-based H/3 indices 
is defined collectively by 74 stars in NGC 6231 
(Shobbrook 1983), NGC 6664 (Schmidt 1982), 
M67 (Nissen et al. 1987) and NGC 6633 (Schmidt 
1976). The standard H/3 indices range in value 
from 2.53 to 2.87, leading to an equation of the 
form /3 = 1.157/3i, with a dispersion of only ±0.015 
mag about the mean line. The mean difference be- 
tween aperture-based H/3 indices from this night 
and the H/3 indices in the profile-fit-based pho- 
tometry set of NGC 6253 was determined with a 
standard error of the mean < 0.001 mag. 

The only Ca frames of NGC 6253 taken on 
a photometric night were obtained on April 22, 
1999. Transformation of the hk indices is tied to 
the observations of field stars found in the cata- 
log of Twarog & Anthony-Twarog (1995). Data 
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from several nights in April 1999 and May 2000 
were combined and analyzed, leading to a trans- 
formation relation of the form hk = 1.05/ifcj. Four 
field star standards observed on April 22, 1999 
established the zero point for the equation cali- 
brating the hk indices specifically for that night, 
with a standard error in the mean for the in- 
tercept of ±0.012. The mean difference between 
the aperture-defined hk indices and the profile-fit- 
bascd hk indices was derived with a standard error 
in the mean value of only ±0.005. Comparable dif- 
ferences between aperture photometry and profile- 
fit-bascd values for other Stromgrcn indices were 
determined with precisions indicated by standard 
errors of the mean for V and b — y of 0.002 mag, 
and 0.003 mag for mi and Ci. 

Calibration oi b — y, mi and Ci indices is gen- 
erally challenging, in part because the parameter 
space of color, abundance, and luminosity class 
must be delicately covered by a wide range of 
standards. Following normal practice, we break 
the calibration equation solutions at a color of 
{b — y)o = 0.42. While past experience has shown 
that cooler dwarfs generally have a slope for the 
b — y transformation relation that is shallower than 
that for the warmer dwarfs and cool giants, no 
cool dwarf members of NGC 6253 are expected to 
be within our survey sample. The transformation 
equation for warmer dwarfs was applied to NGC 
6253 stars with b — y less than 0.62, approximately 
the intrinsic color break of 0.42 with a reddening 
adjustment of 0.20. Transformation equations for 
mi and ci indices for redder stars, presumed to be 
giants, were determined separately based on field 
star standards from the Olsen (1993) catalog. 

Several open clusters observed on May 4 and 5, 
2000 were used to delineate the form of these cal- 
ibration equations. Specific references are NGC 
6231 (Shobbrook 1983), NGC 6633 (Schmidt 
1976), M67 (Nissen et al. 1987) and NGC 3680 
(Nissen 1988). In addition, observations on April 
22, 1999 in IC 4651 were used to calibrate NGC 
6253 indices obtained from aperture photometry 
on that night by comparison to values published 
in Paper I. NGC 6253 was not observed on May 
5; 2000. Field star standards from the catalogs of 
Olsen (1983) and Olsen (1993) were observed on 
both nights in May 2000 and form the basis of the 
intercepts for the calibration equations. 

It became apparent that the conditions for part 



of the night 4 May 2000 must have deteriorated 
from ideally photometric. For example, the stan- 
dard deviations for the calibration equation zero 
point values, based on nine field star standards 
observed on that night, were distressingly large, 
0.04 at best and as high as 0.07 for mi. We con- 
structed an independent check of the calibration 
of Stromgren indices based on the April 22, 1999 
photometry in IC 4651 and several field star stan- 
dards. We were able to estimate corrections to the 
zero points for the transformation equations with 
precisions indicated by the following standard er- 
rors of the mean for the zero points: ±0.002 for V 
and b — y, and ±0.004 for mi and Ci. 

Final photometric values for the stars brighter 
than V ~ 18.3 arc given in Table 1. Pixel coordi- 
nates for the stars are on the (X,Y) system for the 
cluster developed by WEBDA and WEBDA iden- 
tification numbers are provided whenever avail- 
able. The standard errors of the mean as a func- 
tion of V are presented in Fig. 1 for all of the 
indices. The errors are derived from the frame- 
to-frame dispersion for each filter used to con- 
struct the indices rather than the formal errors 
from IRAF. Clearly, the internal errors are excel- 
lent for V < 17.5 for all indices. 

2.3. Comparison to Previous Photometry 

No intermediate-band photometry on any sys- 
tem has been published for NGC 6253. However, 
the V magnitude derived from the y filter can 
be compared to the published broad-band data. 
As detailed in the Appendix, there are systematic 
differences with all the published V photometry, 
though our data arc located near the approximate 
mean of the three photometric surveys. All of the 
broad-band studies exhibit non-negligible offsets 
with respect to each other and, in one case, a 
significant color dependence is found among the 
residuals. We will return to this point in Sec. 5 
when we discuss the cluster match to the theoret- 
ical isochrones. 

3. The Color-Magnitude Diagram: Thin- 
ning the Herd 

The CMD for all stars with V brighter than 
~ 18.3 and at least 2 observations each in b and 
y is presented in Fig. 2. Filled circles are stars 
brighter than V = 17 that have larger than ex- 
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pected variation in V at their magnitude and have 
at least 5 observations in y and multiple observa- 
tions in all filters. The gross features of the clus- 
ter CMD are observable and confirm the results 
found in the broad-band surveys. The turnoff near 
y ~ 15 exhibits a sharp hook to the red below a 
well-defined gap. The subgiant branch is well pop- 
ulated, sloping downward and to the red. The red 
giant clump is easily identifiable near V ~ 12.6, 
with a bright red giant branch that may extend to 
V ~ 11.3 before turning over. The morphological 
resemblance to another definitive old open cluster, 
M67, is striking. 

Unfortunately, due to its galactic location, con- 
tamination of the CMD by non-members is a dom- 
inant problem for any analysis. Below V = 16.5, 
the cluster main sequence merges into the amor- 
phous band defined by background field stars, 
while field subgiants and giants dominate for b — y 
redder than 0.7 and V > 15. The problem is simi- 
lar in nature to that faced in Paper I for IC 4651, 
with the complication that the cluster turnoff lies 
much closer in color and magnitude to the field 
star distribution due to the greater distance of 
NGC 6253. 

As discussed in Paper I, the focus of our analy- 
sis of this cluster, and all clusters in the program, 
is the precise determination of the primary clus- 
ter characteristics of reddening and metallicity fol- 
lowed, when possible, by distance and age. Pho- 
tometric derivation of these parameters requires 
color indices of high precision for single-star mem- 
bers of the cluster. If we have achieved the photo- 
metric accuracy we require, the dominant issue be- 
comes the isolation of highly probable members of 
the cluster from the field distribution. Procedures 
that identify and eliminate probable non-members 
arc acceptable even if they remove some cluster 
members. Keeping non-members out is more im- 
portant than keeping members in since complete- 
ness of the sample is not a issue. In the sections 
that follow, we will refine the sample in Fig. 2 to 
reach a point where the probable cluster members 
alone define the features of the CMD. 

3.1. Thinning the Herd: The Radial Dis- 
tribution 

In the absence of kinematic information to es- 
tablish cluster membership, the next best option 
is to use the radial distribution of cluster stars 



to select the area with the highest proportion of 
members. Previous photometric surveys have es- 
tablished estimates of the cluster diameter, begin- 
ning with an estimate of 5' from Lynga (1987). 
The limiting field size of CCD's has been a factor 
in subsequent photometric studies. BR studied 
two overlapping fields that were 6.8' x 6.8', one 
centered on the cluster and one north of center. 
SA imaged a cluster field 2.1' x 3.3' along with a 
comparison field 10' to 15' away from the cluster 
center. PI observed a field approximately 4' x 4'. 
With the exception of BR, the fields studied have 
been smaller than the proposed cluster diameter 
and, even in the case of BR, the counts beyond 
the cluster core cover a small portion of the area 
around the cluster. The location of the cluster 
center has been selected by a simple visual fit in 
the case of BR, while PI and SA fixed the cen- 
ter by adjusting the coordinates to maximize the 
cluster density profile. The locations found by BR 
and SA are similar, but differ from that of PI by 
almost 1'. 

With a larger field size, we hoped to have an 
easier time establishing the profile and center of 
the cluster. Coordinates for stars brighter than 
limits varying from Yum = 15.0 to 18.0 were used 
to determine the cluster center which would op- 
timize the concentration of the cluster distribu- 
tion within the smallest radius. The cluster center 
we found is similar to that found by BR and SA, 
though the uncertainty in the exact location seems 
to be larger than claimed by SA, due to the small 
number of bright stars defining the center. The 
cluster blends into the background at a diameter 
of 7', larger than foimd by Lynga (1987). 

Due to the modest portion of the CCD field oc- 
cupied by the cluster, the dominance of the CMD 
by non-members is easily understood. The clus- 
ter excess above background is even smaller for 
V > 17.5, showing an increasing deficiency of 
lower mass stars compared to expectations from 
normal luminosity functions. To maximize the 
percentage of cluster members, our first cut will 
remove all stars located more than 100" away from 
the cluster center. Simultaneously, we will identify 
all stars that have standard errors in the mean 
for b — y greater than 0.010 for later elimination. 
The resulting CMD is presented in Fig. 3, where 
crosses identify stars that fail to meet the error 
cut. 
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The primary features of the cluster CMD are 
significantly enhanced, though confusion caused 
by field star contamination remains an issue be- 
low V = 17. Beyond the features noted earlier, 
one can now identify a binary sequence located 
about 0.7 mag above the expected position of the 
main sequence. The exact location of the main 
sequence is still a source of controversy due to the 
band of field stars that extends from b — y = 0.55 
the region redward of the subgiant branch. 

3.2. Thinning the Herd: CMD Deviants 

Given the high precision of the b — y indices, 
if the cluster sample dominated that of the field 
stars, one could define the mean relation for the 
main sequence in the CMD and exclude all stars 
which deviated by more than three sigma from the 
mean, eliminating both field stars and probable 
binaries. In this situation, where the field stars 
contribute as much, if not more, to the region of 
the main sequence, definition of the cluster main 
sequence can become an exercise in selection bias. 

Fortunately, the availability of multiple indices 
makes the selection process much less subjective. 
To demonstrate the point, we make use of the fil- 
ter pair with the largest baseline, u — y. For our 
next cut, we include stars with V brighter than 
17.25 and bluer than b — y = 0.7, at least 2 ob- 
servations each in y, h, v, and u and a standard 
error of the mean in Ci below 0.025 mag. The 
last constraint guarantees that the error in u — y 
is below 0.020. The (y,u - y) CMD is shown in 
Fig. 4. The increased temperature sensitivity of 
w — y is obvious and allows easy identification of 
stars that are too blue for their apparent magni- 
tude (field stars), plotted as open triangles in the 
figure, and stars that are too red (binaries or field 
stars), plotted as crosses. The success of the pro- 
cess is measurable in Fig. 5, where the traditional 
CMD for the stars is presented. With one excep- 
tion, the stars that lie within the binary sequence 
in the standard CMD have been identified as such 
from the u — y index. Elimination of the stars 
that scatter blueward in the u — y figure would 
remove all the scatter to the blue in the standard 
CMD. Though it is probable that elimination of 
these deviants from the analysis will sweep away a 
few true cluster members, the critical point is that 
the remaining sample contains a high percentage 
of single stars whose structure and evolution are 



typical of the cluster as a whole. 

4. Fundamental Properties: Reddening 
cind Metallicity 

The select sample of probable members (filled 
circles) in Fig. 5 includes 99 stars. Our final re- 
striction is to impose a color cut on the turnoff. 
Only stars in the color range from b — y = 0.53 to 
0.60 will be used to determine the reddening and 
metallicity. The blue limit removes the blue strag- 
glers from the analysis. Since these stars have, by 
definition, anomalous evolution, conclusions based 
upon their photometry are always open to ques- 
tion. The red cutoff removes the evolved stars of 
the subgiant branch and keeps the sample within a 
color range where the intrinsic calibrations of the 
uvby system can be expected to produce reliable 
parameters. The final sample includes 71 stars, 
none of which have errors in H/3 above ±0.020; 
only 2 stars have photometric uncertainty in H/3 
above ±0.015. 

4.1. Reddening 

As discussed in Paper I, derivation of the red- 
dening from intermediate-band photometry is a 
straightforward, iterative process given reliable es- 
timates of H/3 for each star. Briefly, the intrin- 
sic 6 — y color is derived from the reddening and 
metallicity-indepcndent H/3 index. Since b—y is af- 
fected by both metallicity and evolutionary state, 
the intrinsic value is then adjusted to correct for 
the differences between the star and the unevolved 
stars used to define the standard relation linking 
b — y and H/3. Because the indices estimating the 
metallicity and luminosity are influenced by the 
reddening, one must iterate the procedure, cor- 
recting each stellar index for the derived redden- 
ing, calculating a new set of intrinsic indices, lu- 
minosities, and metallicities, then re-deriving a re- 
vised reddening. The procedure converges quickly. 
For cluster stars, all members are assumed to have 
the same reddening and metallicity, so one may 
calculate the reddening under different fixed as- 
sumptions for metallicity, and the metallicity un- 
der different fixed assumptions for the reddening. 
Only one combination of reddening and metallic- 
ity will produce mutually consistent calculations. 

Since both procedures give effectively identical 
results, the primary decision is the choice of the 
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standard relation for H/? versus b — y and the ad- 
justments required to correct for metahicity and 
evolutionary state. The two most commonly used 
relations are those of Olscn (1988) and Schuster 
& Nissen (1989). As found in Paper I for IC 
4651, both produce very similar but not identi- 
cal results. We have processed the indices for the 
71 stars through both relations and find E{b — y) 
= 0.192 ± 0.017 (s.d.) with Olsen (1988) and 
E{b -y) = 0.188 ± 0.015 (s.d.) with Schuster & 
Nissen (1989). As a compromise, we will take 
the weighted average of the two and use E{b — y) 
= 0.190 ± 0.002 (s.e.m.) or E{B -V) = 0.260± 
0.003 (s.e.m.) in the analysis that follows. 

It seems useful to discuss the reliability of the 
reddening estimate before deriving the metallic- 
ity in the next section. Though we have used the 
turnoff stars alone in the calculation, we can check 
the consistency of the final value using all the stars 
in the members sample, i.e., by including the blue 
stragglers. Fig. 6 shows a plot of the standard re- 
lation between and h — y, reddened by E{b — y) 
= 0.190 and superposed upon the observed data 
for NGC 6253. With the exception of 3 blue strag- 
glers (indicated by filled symbols) that lie system- 
atically above the standard relation, the adjusted 
curve falls nicely through the mean of the hot- 
ter points defined by the A through early F stars. 
Note that the curve does not pass through points 
defined by the stars at the turnoff. The reason 
is that the metallicity effect on the cooler stars is 
significant and no metallicity correction has been 
applied to the cooler end of the standard relation. 
The systematic difference relative to an approxi- 
mately Hyades metallicity is the first indication of 
the exceptional metallicity of the cluster. 

As for the three stragglers located above the 
mean relation, the significantly lower reddening es- 
timate implied for these stars is explained if they 
are foreground stars rather than members of the 
cluster. This is confirmed by ci indices for these 
stars which are systematically smaller by approx- 
imately 0.1 mag than the expected values at their 
b — y based upon the rest of the stragglers. If we 
apply the cluster reddening to these stars, the cor- 
rected indices would generate physically implausi- 
ble combinations of (6 — y)o and cq. Moreover, 
use of the individual reddening values for these 
stars rather than the adopted cluster value pro- 
duces metallicity estimates close to those typical 



of stars in the solar neighborhood, rather than the 
exceptionally high value found for the cluster. 

How does the derived reddening value compare 
with other estimates? At the galactic coordinates 
of the cluster, the reddening map of Schlegel et 
al. (1998) based upon infrared emission by dust 
gives E{B — V) = 0.35. Since the cluster is lo- 
cated more than 150 pc away from the galactic 
plane and the reddening map includes all redden- 
ing along the line of site at and beyond the cluster, 
the difference between the limit and the observed 
value is understandable. 

The derived value falls nicely within the range 
of E{B -V) = 0.23 to 0.32 found by BR using 
comparisons with theoretical cluster CMDs but 
the agreement may be fortuitous if the problems 
with the broad-band photometry discussed in the 
Appendix are tied to the data of BR. The revised 
values of B — V from BR transformed to the sys- 
tems of SA and PI expand the color difi^erential 
between the giant branch and the turnoff, inde- 
pendent of issues tied to the zero-points, leading 
to a revision of what age and metallicity combina- 
tions supply the optimum fit to the isochrones. 

Both PI and SA use multicolor CMDs matched 
to standard sequences (PI) and theoretical isochrones 
(SA) to derive E{B -V) = 0.20 and E{V - I) of 
0.25 to 0.27. PI quote uncertainties of ±0.05 mag 
based upon the internal scatter of the main se- 
quence alone, while SA provide no error estimate. 
The fact that both investigations produce simi- 
lar estimates of E{B — V) and E{V — I) may be 
purely fortuitous. PI matches the cluster to a. BV 
sequence tied to stars typical of the solar neigh- 
borhood, i.e., metallicity less than or equal to so- 
lar, while SA fits isochrones with Z = 0.05, more 
that double solar metallicity. The VI diagrams 
should be relatively insensitive to metallicity ef- 
fects, but comparison of SA and PI demonstrates 
a zero-point discrepancy in their {V — I) colors of 
0.14 mag. In summary, while realistic appraisal of 
the previous cluster reddening analyses indicates 
an estimate of E{B — V)= 0.2 ±0.1, this result 
doesn't offer much value in testing the reliability 
of the intermediate-band result. 

4.2. Metallicity from mi 

Given the reddening of E{b — y) = 0.190, 
the derivation of [Fe/H] from the mi index is 
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a straightforward procedure. The mi index for 
a star is compared to the standard relation at 
the same color and the difference between them, 
adjusted for possible evolutionary effects, is a mea- 
sure of the relative metallicity. Though the com- 
parison of mi is most commonly done using b — y 
as the reference color because it is simpler to ob- 
serve, the preferred reference index is H/? due to 
its insensitivity to both reddening and metallicity. 
Changing the metallicity of a star will shift its 
position in the mi - (b — y) diagram diagonally, 
while moving it solely in the vertical direction in 
m-i - H/3. Moreover, reddening errors do not lead 
to correlated errors in both mi and H/3. 

The primary weakness of metallicity determi- 
nation with intermediate-band filters is the sensi- 
tivity of [Fe/H] to small changes in mi ; the typical 
slope of the [Fe/H], Sirii relation is 11. Even with 
highly reliable photometry, e.g., mi accurate to ± 
0.015 for a faint star, the uncertainty in [Fe/H] for 
an individual star is ±0.17 dex from the scatter 
in mi alone. When potential photometric scat- 
ter in H/3 and ci are included, errors at the level 
of ±0.2 dex are common. As noted in previous 
papers in this series, the success of the technique 
depends upon both high internal accuracy and a 
large enough sample to bring the standard error of 
the mean for a cluster down to statistically useful 
levels, i.e., below 0.10 dex. 

Because of the size of the sample, we can mini- 
mize the impact of individual points such as bina- 
ries and/or the remaining non- members, though 
they will clearly add to the dispersion. After cor- 
recting each star for the effect of E(b — y) = 0.190, 
the m,i, ci, and H/3 values for the 71 core cluster 
members were inn through the metallicity calibra- 
tion using the relation derived by Olsen (1988) 
and also applied to the metal-rich open cluster IC 
4651 in Paper I. Before discussing the results, it 
should be noted that the range in Hfi among the 
turnoff sample is small. Since all the stars have 
approximately the same color, variations in mi 
should be primarily a product of the evolutionary 
state, normally measured via (5ci, the difference 
between the observed ci and the standard value at 
the same color. As demonstrated in Nissen et al. 
(1987) from photoelectric observations of M67, for 
the cooler F stars, evolution reduces the observed 
value of mi at a given color, making the more 
evolved stars appear more metal-deficient. Since 



this effect occurs exactly within the color range of 
our sample, we have adjusted the mi indices of 
the more evolved stars upward by 0.3Aci, where 
Aci is the mean ci for a star minus the average 
ci value of the stars fainter than V = 16.5. We 
emphasize that while the brighter stars do, in the 
mean, have smaller mi and larger ci, the correc- 
tion implied by the above procedure is small and 
does not completely eliminate the correlation be- 
tween m,i and V . 

From 71 stars, the average 5mi is -0.071 ± 
0.002 (s.e.m.), which translates into [Fe/H] = 
±0.90 ± 0.03 on a scale where the Hyades has 
[Fe/H] = ±0.12. The 5mi index and its implied 
metallicity are unmatched by any star or cluster 
observed to date within the Milky Way. Though 
the extreme nature of the estimate casts doubt 
on its validity, before evaluating what this means, 
we turn to another metallicity indicator, hk, as a 
means of constraining our potential solutions. 

4.3. Metallicity from hk 

The hk index is based upon the addition of 
the Ca filter to the traditional Stromgren fil- 
ter set, where the Ca filter is designed to mea- 
sure the bandpass which includes the H and K 
lines of Ca 11. The design and development of 
the Caby system have been laid out in a se- 
ries of papers discussing the primary standards 
(Anthony-Twarog et al. 1991b), an extensive cat- 
alog of field star observations (Twarog & Anthony- 
Twarog 1995), and calibrations for both red giants 
(Anthony-Twarog & Twarog 1998) and metal- 
deficient dwarfs (Anthony-Twarog et al. 2000). 
Though the system was optimally designed to 
work on metal-poor stars and most of its ap- 
plications have focused on these stars (Anthony- 
Twarog et al. 1995; Baird 1996), early indications 
that the system retained its metallicity sensitivity 
for metal-rich F dwarfs have been confirmed by 
observation of the Hyades and analysis of nearby 
field stars (Anthony-Twarog et al. 2002). 

What makes the hk index, defined as [Ca — 6) — 
(b — y), so useful even at the metal- rich end of the 
scale is that it has half the sensitivity of mi to 
reddening and approximately twice the sensitivity 
to metallicity changes, i.e., A[Fe/H]/A/ifc =5.6. 
The metallicity calibration for F stars derived in 
Anthony-Twarog et al. (2002) used 6hk defined 
relative to 6 — y as the temperature index. To 
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minimize the impact of reddening on metallicity, 
this calibration has been redone using H/3 as the 
primary temperature index. The Hyades hk,{b—y) 
relation has been converted to hk, H/3 through the 
standard (b — y), H/3 relation discussed in Sec. 4.1 
and illustrated in Fig. 6. Prom 50 stars compiled 
in Anthony- T war og et al. (2002), one finds 

[Fe/H]= -3.516hk{f3) + 0.12 

with a dispersion of only ±0.09 dcx about the 
mean relation. Though the derived zero-point of 
the relation was found to be +0.07, it has been 
adjusted to guarantee a Hyades value of +0.12, 
the same zero-point used for the mi calibration. 

We have derived the cluster metallicity using 
the hk indices for the same 71 stars analyzed above 
using both hk relative to the b — y relation and hk 
relative to H/3. In order to be consistent with our 
approach to mi, we have included an evolution- 
ary correction to hk of the same manner and size 
adopted for ^toi . The necessity for such a correc- 
tion is confirmed by a plot of the mean hk versus V 
where one finds that the more evolved stars have 
systematically lower values of hk, though the size 
of the trend is comparable to that found for mi. 
The results relative to (6 — y) and H/3 are [Fe/H] 
= +0.84 ± 0.14 (s.d.) and +0.68 ± 0.14 (s.d.), 
respectively. The first point to note is the signif- 
icantly smaller size of the dispersion for hk-h&sed 
[Fe/H] in comparison with mi, despite that fact 
that there are fewer Ca frames than v frames used 
in the construction of the indices and the errors in 
hk are comparable to or larger than those for mi . 
This is a direct benefit of the smaller slopes in the 
calibration curves. 

Second, the [Fe/H] from the (6 — 2/)-based rela- 
tion is larger by a significant amount, though simi- 
lar to the mi abundance. The difference is a reflec- 
tion of the high sensitivity of the [Fe/H] to small 
changes mb — y, particularly through the adopted 
reddening. If we lower the reddening by 0.02, i.e., 
increase the intrinsic b — y hy the same amount, 
the (6-j/)-based [Fe/H] declines by 0.22 dcx. The 
same reddening adjustment only alters the H/3- 
based [Fe/H] by 0.01 dex. Thus, demanding that 
both approaches give identical results would re- 
quire a shift of about -0.015 in E{b — y), not an 
implausible change given the photometric uncer- 
tainty. A shift of this size only alters the wi-based 
[Fe/H] by 0.07 dex, producing an average for hk 



and mi of [Fe/H] = +0.75. 

Third, though the analyses of mi and hk lead 
to [Fe/H] ranging from +0.65 to +0.9, an annoy- 
ing range given the photometric accuracy, there is 
complete qualitative agreement in the sense that 
both indices imply that this cluster is more metal- 
rich than any cluster or object ever studied on 
either system. Assuming that the extreme nature 
of the indices is not the product of photometric 
error, the unique [Fe/H] estimates derived by ex- 
trapolation of the standard calibrations brings the 
validity of their absolute numerical value into obvi- 
ous question. In short, NGC 6253 resides in unex- 
plored photometric territory and the initial foray 
into this strange land leads to the conclusion that 
it cannot be regarded as minor variant on what we 
are familiar with locally. Since both photometric 
systems have been tested for stars up to [Fe/H] 
= +0.3, it seems reasonable to assume that the 
cluster abundance is, at least, this high. As a sim- 
ple compromise, we will assume that [Fe/H] lies 
somewhere between +0.4 and +0.6. 

Before discussing the age and distance deter- 
mination, the effects of errors in the photometry 
should be investigated. It has already been noted 
above that a small shift in the reddening would 
help to bring the photometric results into better 
agreement, though still at a high value of [Fe/H]. 
Small errors in E{b — y) can result from indepen- 
dent errors in H/3, b — y or both, or from a failure 
of the intrinsic color calibration at high [Fe/H]. 
Likewise, one could make the claim that the high 
metallicity is a product of errors in the zero-points 
for the mi and hk transformations. 

The weakness in most explanations that at- 
tempt to resolve the question by citing an error in 
an individual index or parameter is that they only 
resolve a portion of the problem. For example, 
a zero-point error in mi has little impact on hk. 
Adjusting the reddening due to a claim that the 
intrinsic color relation underestimates the metal- 
licity effect on b—y would significantly alter metal- 
licity indices tied to b—y, but have little impact on 
those based upon H/3. The only exception to this 
rule would be an error in the zero-point for the 
H/3 transformation. If we adjust the H/3 indices 
by -0.030, the cluster reddening declines by the 
same amount. The combined effect on the metal- 
licity is to lower [Fe/H] from mi (H/3), /ifc(H/3), 
and hk{b-y) to +0.47, +0.33, and +0.47, respec- 
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tively. The agreement among the indices is better 
and at a somewhat lower metaUicity, but it still 
places NGC 6253 above the current limit of the 
metallicity scale. We also emphasize that every 
test we have made to check the zero-point of the 
H/3 photometry, night-to-night comparisons, stan- 
dard star observations that bracket the cluster on 
a given night, indicate that if an offset is required, 
it must be small and it must make H/3 larger, not 
smaller. 

5. Fundamental Parameters: Distance 
cmd Age 

5.1. The Broad-Band CMD 

Given the reddening and metallicity, the tradi- 
tional approach to deriving the distance and age is 
to compare the cluster CMD to a reliably defined 
set of theoretical isochrones. By reliably defined 
we mean color transformations and bolometric 
corrections between the theoretical and the obser- 
vational plane that reproduce the colors and abso- 
lute magnitudes of nearby stars with known tem- 
peratures and abundances. From an age stand- 
point, the critical test is whether or not a star 
of solar mass, composition, and age resembles the 
sun. For open clusters of solar and subsolar abun- 
dance, the impact of this issue on cluster ages 
and distances has been emphasized on a variety 
of occasions (Twarog & Anthony- T war og 1989; 
Twarog ct al. 1993; Daniel et al. 1994; Twarog 
et al. 1995, 1999). 

For NGC 6253, two problems immediately 
arise. First, isochrones are invariably created for 
broad-band systems and a check of the most recent 
publications shows that the theoretical isochrones 
are available on the U BVRI system, among oth- 
ers, but not uvby. Second, though isochrones have 
been generated with a mass fraction of metals up 
to five times the solar value, the questions of the 
related helium fraction and whether or not the 
metals should be scaled-solar or enhanced in the 
a-elements remain open. Moreover, the lack of 
a significant population of such stars with reli- 
ably determined fundamental parameters makes 
the color transformations an exercise in faith in 
the validity of extrapolating synthetic colors from 
model atmosphere codes. 

Fortunately, the first problem is easily solved. 
The b — y index has long been known to corre- 



late well with B — V at sl given [Fe/H] with little 
dependence on evolutionary state. Though one 
could use the published BV data to construct a 
CMD of cluster members, there are zero-point dif- 
ferences among the broad-band surveys, the pho- 
tometric accuracy of the surveys is not uniform, 
and the surveys only include stars within the core 
of the cluster. Instead, we chose to transform our 
b — y data to B — V. The construction of the 
mean broad-band photometric indices through the 
merger of the published data is discussed in the 
Appendix; for reasons noted there, we have cho- 
sen the B — V zero-point and scale of SA as the 
base for the merger. Using B — V and 6 — y for 
all the stars brighter than V = 16.5 common to 
the merged BV sample and our by data, linear re- 
lations were derived between the two indices. To 
optimize the color transformation, the sample was 
subdivided into three color regimes, with the tran- 
sition points at b — y = 0.47 and 0.71. Exclud- 
ing 4 stars with residuals about the mean relation 
greater than 0.10 mag, the transformation curves 
over the three color regimes are: 

b-y< 0.47, N=30 
B-V = 1.492(±0.038)*(6 -y) + 0.009(±0.016) 

0.47 <ib-y)< 0.71, N=264 
B-V = 1.717(±0.026)*(6 -y)- 0.096(±0.015) 

{b-y)> 0.71, N=85 
B-V = 1.464(±0.029)*(6 -y) + 0.084(±0.024) 

The standard deviations of the residuals about 
the mean relation for the three regions are 0.017, 
0.020, and 0.017 mag, respectively, consistent with 
a typical uncertainty of 0.015 mag in B — V and 
less than 0.010 mag in 6 — y. A plot of the data 
and the mean relations is shown in Fig. 7. 

To define the cluster CMD, we restrict our sam- 
ple to stars within the core for 13.0 < < 17.25 
with a standard error in the mean for b ~ y be- 
low 0.010 (see Fig. 3). The only additional ex- 
clusion will be the stars identified in Fig. 4 by 
triangles and tagged as nonmembers based upon 
a. u — y index that places them blueward of the 
cluster main sequence. For V < 13.0, all stars ob- 
served within the entire field will be used in an 
effort to strengthen the definition of the red giant 
clump and the giant branch. 
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5.2. The Isochrones 

The two primary sets of theoretical isochrones 
available in an array of broad-band systems, 
metallicities, and ages are those of Girardi et al. 
(2002) (hereinafter referred to as PAD) and Leje- 
une & Schaerer (2001) (hereinafter referred to 
as GEN). Based upon the many comparisons, in- 
cluding our own (Anthony-Twarog et al. 1994; 
Daniel et al. 1994; Anthony-Twarog et al. 1991a; 
Twarog et al. 1995), between open clusters and 
past and present generations of isochrones, we will 
only make use of isochrones that include convec- 
tive overshoot mixing. On a scale where solar 
metalhcity is Z 0.019 and Y = 0.273, PAD 
supply isochrones with (Y,Z) = (0.32,0.040) and 
(0.39, 0.070) with either scaled-solar abundances 
or a-enhanced abundances. In contrast, the GEN 
set assumes solar (Y,Z) = (0.30,0.02) and has 
isochrones available for (0.34,0.04) and (0.48,0.10). 
While the latter metallicity is closest to the clus- 
ter abundance implied by the photometry, the 
combination with the high helium abundance pro- 
duced a CMD morphology that was so discrepant 
in comparison with the cluster that inclusion of 
the set added nothing to the discussion. While in 
the past we have attempted to test the zero-points 
of the isochrones by comparison to the sun and/or 
field stars (Twarog et al. 1999), there is no way of 
knowing if a difference found between theory and 
observation at solar metallicity transfers directly 
to isochrones at [Fe/H] = -1-0.3 to -1-0.6. Because 
we are fixing the reddening at E{B — V) = 0.26, 
differences among the derived ages and distances 
are solely a product of the differences among the 
models, both in assmned composition and in the 
transformations between the theoretical and the 
observational plane. 

For the first isochrone fit, we use the GEN mod- 
els for (Y,Z) = (0.34, 0.04), equivalent to [Fe/H] 
= -1-0.3, the low end of the metallicity range po- 
tentially applicable to NGC 6253. The isochrones 
in Fig. 8 have been shifted by -1-0.26 in {B - V) 
and by 11.60 in V, the apparent distance modu- 
lus. Filled circles and crosses are the same core 
members and probable binaries identified in Fig. 
4. Open triangles are the additional data from 
the core of the chister outside the c;olor range of 
Fig. 4 and squares are the points outside the 
core brighter than V = 13.0. The superposed 
isochrones have ages of 3.2, 3.55, and 4.0 Gyrs. 



The match to the isochrones at the turnoff and 
the giant branch is very good, though the curva- 
ture of the main sequence turnoff is slightly larger 
than predicted by the models. Because of the in- 
clusion of the brighter stars over the entire field, 
the giant branch can now be mapped to the mod- 
els well above the level of the clump at V = 12.7. 
The overall fit of the data implies an age between 
3.4 and 4.0 Gyr, so we will adopt 3.7 ± 0.3 as the 
value for these models. 

The next set of models will be those of PAD 
with scaled-solar abundances and (Y,Z) = (0.32, 
0.04) or [Fe/H] = -1-0.32. Despite the similarities 
in [Fe/H] to the isochrones of Fig. 8, a match to 
the turnoff and main sequence requires an appar- 
ent modulus of 12.30 in Fig. 9. The ages of the 
superposed isochrones are 2.0, 2.5, and 3.2 Gyrs. 
The overall quality of the fit is significantly re- 
duced. The turnoff exhibits the same difference 
in curvature found in Fig. 8, but the subgiant 
branch lies between the two older tracks while the 
color of the turnoff is between the two younger. 
More important, the color of the isochrone giant 
branches is too red by about 0.2 mag inB — V and 
the red giant clump is too faint compared to the 
observations by 0.6 mag. If one reduces the red- 
dening to E{B — V)= 0.20, the apparent modulus 
is lowered to 11.95, producing a consistent match 
between the subgiant branch and the turnoff color, 
but leaving the clump too faint by 0.25 mag and 
the giant branch of the isochrones too red by about 
0.08 mag. 

The next alternative is the PAD models with 
(Y,Z) = (0.32,0.04) but enhanced-a elements. The 
results are shown in Fig. 10 for ages of 2.0, 2.5, and 
3.2 Gyrs and an apparent modulus of 12.05. In 
contrast with the previous isochrones, the match 
from the main sequence through the base of the gi- 
ant branch is effectively perfect in magnitude and 
color. The curvature of the turnoff is exactly as 
predicted. The giant branch of the isochrones is 
still too red by between 0.03 and 0.06 mag, but the 
clump is much closer in luminosity to the observed 
data. The age of the cluster is easily constrained 
to be 2.5 ± 0.1 Gyr. 

The last match shown in Fig. 11 is that 
for the PAD models, a-enhanced, with (Y,Z) = 
(0.39,0.07) or [Fe/H] = -^0.57. The adopted ap- 
parent modulus is 12.15 and the isochrones have 
ages of 1.0, 1.4, and 2.0 Gyr. As in Fig. 9, 
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the quality of the fit is disappointing, particularly 
since this represents an overall metallicity closest 
to that defined by the photometric indices. The 
observations of the turnoff and subgiant branch 
cannot be reconciled with a single age and the dis- 
crepancy between the locations of the giant branch 
is large, with the usual pattern of the models be- 
ing too red and the theoretical clump too faint. As 
with Fig. 9 one could slightly improve the turnoff 
match by lowering the reddening and shifting an 
older isochrone into alignment with the cluster, 
but the subgiant branch would still be found at a 
magnitude level inconsistent with the age defined 
by the color of the turnoff. Use of the scaled-solar 
isochrones at Z = 0.07 produces an increase in 
the apparent modulus to 12.55 and an even worse 
match to the shape of the turnoff and the position 
of the giant branch. 

Given the large range of parameters derived us- 
ing the various isochrones, can we reach any plau- 
sible conclusions regarding the cluster age and dis- 
tance? The first point to note is that, with the 
exception of the Z = 0.07 isochrones, the models 
used were of lower Z than implied by strict accep- 
tance of the photometric metallicities. Increas- 
ing Z for a cluster at a fixed reddening leads to 
a younger age and a larger distance modulus, as 
confirmed above. If the age and reddening were 
known, one could derive the metallicity by se- 
lecting the isochrones that optimally matched the 
CMD, a simplified version of the approach taken 
by BR. 

As discussed earlier, the morphology of the 
CMD, in particular the existence of a main se- 
quen('(^ hook and gap, constrains the cluster to be 
younger than about 5 Gyr. A more quantitative 
answer can be attained by using the morphologi- 
cal age ratio (MAR) as discussed and revised by 
Anthony- Twarog & Twarog (1985) and Twarog 
& Anthony- Twarog (1989). Based upon the ra- 
tio of the difference in V between the red giant 
clump and the brightest point at the turnoff and 
the difference in i? — V between the giant branch 
at the level of the clump and the vertical turnoff, it 
has the advantages of being more sensitive to age 
changes than either difference alone, reddening- 
independent, and less sensitive to metallicity than 
either difference alone. From the CMD of NGC 
6253, the MAR value is found to be 2.9 (2.8 if one 
adjusts for the possible difference in reddening for 



the giants compared to the turnoff stars (Twarog 
et al. 1997)), slightly smaller than M67 at 3.25 
but larger than NGC 6819 at approximately 2.0. 
If we adopt ages of 4.0 Gyr (Dinescu et al. 1995) 
and 2.5 Gyr (Rosvick & Vandenberg 1998) for 
these clusters, the approximate age of the NGC 
6253 is 3.5 Gyr, in very good agreement with the 
values typically derived by BR from a more in- 
volved match to a variety of isochrones. Keeping 
in mind the caveat that we have no knowledge of 
how reliable an index such as MAR is for [Fe/H] 
of -1-0.4 or higher, the GEN isochrones are con- 
sistent with this age for the cluster if [Fe/H] is 
increased to -1-0.4. At this metallicity, the appar- 
ent modulus becomes 11.7, leading to My for the 
clump of +1.0, consistent with expectation for an 
older cluster of high metallicity as empirically de- 
termined by Twarog et al. (1997). 

A unfortunate issue that we cannot resolve 
is the disagreement between the scaled-solar 
isochrones of GEN and PAD. Though the GEN 
models use Y higher by 0.02, the differences in age 
(~ 1 Gyr) and apparent modulus (0.7 mag) seem 
excessive given the similar values of Z. The lower 
age defined by the turnoff color is consistent with 
the discrepancy between the colors of the observed 
and computed giant branches; the observed CMD 
morphology requires an older age than allowed by 
the absolute colors. This explains why artificially 
lowering the reddening produces a better match 
to an older isochrone and why increasing Z to 0.07 
makes the problem worse. Following the pattern, 
the predicted My values of the clump of -1-0.4 
and -1-0.15 for solar-scaled Z = 0.04 and 0.07, re- 
spectively, are inconsistent with a very metal-rich 
cluster at these ages. 

Despite the problems, the PAD isochrones may 
provide a solution to the entire issue. By far, 
the best match to the shape of the observed 
turnoff and subgiant branch was supplied by the 
a-enhanced isochrones with Z — 0.04. The im- 
provement over any solar-scaled set may indicate 
that the true [Fe/H] for NGC 6253 is closer to 
-f 0.4, but the photometric indices, mi and hk, give 
anomalous abundances because the energy distri- 
butions of the stars are distorted by the enhanced 
a elements. The morphology of the CMD and the 
My of the clump cannot be compared reliably to 
nearby clusters of the same age because these clus- 
ters have predominantly scaled-solar distributions. 
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The My of the clump with the optimal match of 
Fig. 10 is +0.65, within the observed range for 
nearby clusters of the same age and close to the 
value predicted by the isochrones. 

6. Summary and Conclusions 

The open cluster, NGC 6253, has unusual sig- 
nificance in two ways. It populates a rather rare 
class of open clusters, old and significantly closer 
to the center of the Galaxy than the sun. The 
statistical studies of Janes & Phelps (1994) and 
Hufnagel & Smith (1994) have demonstrated that, 
observational difficulties aside, the number of open 
clusters observed within 7 kpc of the galactic cen- 
ter is so small that the probability for survival be- 
yond 1 Gyr is low; NGC 6253 falls just outside 
this zone. Second, only the old open cluster NGC 
6791 has been identified as having a metallicity 
significantly higher than the Hyades (Chaboyer et 
al. 1999), and even this claim has been questioned 
(Taylor 2001). Past attempts to quantify the key 
parameters for NGC 6253 have depended primar- 
ily on consistency fits to multicolor, broad-band 
relations and synthetic CMD's tied to theoretical 
isochrones. Though there is a consensus linked to 
the CMD morphology that the cluster must be 2 to 
4 Gyr old and more metal-rich than the sun, nar- 
rowing the range has been difficult. Though the 
source of the disagreement has originated mostly 
with the different predictions from the various 
isochrones, a problem repeated in Sec. 5, the limi- 
tations imposed by the differences among the pho- 
tometric surveys must be accounted for. As noted 
in the Appendix, the absolute answer one derives 
for the parameters would depend upon which zero- 
point and which color slope one chooses as the 
correct system for defining the broad-band CMD 
and color-color plots. 

We reiterate one of the chief advantages of the 
intermediate-band approach because it allows to 
fix the reddening through H/3 with only a weak 
dependence on the cluster metallicity. Removing 
this as a free parameter in any match to theoret- 
ical isochrones significantly restricts the allowed 
combinations of metallicity, distance, and age. As 
has now been demonstrated using IC 4651 (Paper 
I), NGC 6397 (Paper II), and NGC 6253, it is pos- 
sible to derive the distance with an internal error 
of only ±0.1 mag, with internal errors of the mean 



for E{b-y) and [Fe/H] of ±0.002 and ±0.02, even 
though we cannot guarantee that all non- members 
and binaries have been eliminated. 

Given this precision, how does one explain the 
extraordinary metallicities implied by both the mi 
and hk indices? The flip side to the internal preci- 
sion that is achievable is that the primary source 
of uncertainty in deriving photometric-based pa- 
rameters from a large statistical sample falls upon 
the transformations of the instrumental indices to 
the standard system. The challenge of doing this 
well, even for an established broad-band system 
like UBV, is obvious from the comparisons of the 
data on NGC 6253 as given in the Appendix. How- 
ever, as noted in Sec. 4, to simultaneously reduce 
the cluster metallicity by claiming an error in the 
photometry requires either a conspiracy between 
mi and hk or an error in the common link for the 
reddening and the both metallicity indicators, H/3. 
Unfortunately, an error of over 0.03 mag would be 
required to have a significant impact on the metal- 
licity and, even for this change, would produce 
[Fe/H] well above the Hyades. 

A probable solution to the question is that the 
indices are, in fact, close to the truth but their in- 
terpretation is based upon a questionable extrapo- 
lation of the standard relations for stars with more 
"normal" abundances, both in an absolute and a 
relative sense. This distortion could apply to both 
the metallicity and the reddening determination, 
though the size and direction of the effect on the 
latter variable are unknown. The existence oluvhy 
index distortions due to abundance anomalies has 
been known for over 20 years with CN variations in 
CH and Ba stars producing anticorrelated changes 
in mi and ci being a prime example (see, e.g.. 
Bond 1980; Luck & Bond 1982). The correlated 
impact of abundance anomalies on m,i, ci, and 
hk has been documented for field stars (Anthony- 
Twarog et al. 1991b, 1992) and globular clusters 
like M22 (Anthony- Twarog et al. 1995) and w 
Cen (Rey et al. 2000), but the objects under 
discussion have been metal-poor giants contained 
within systems that exhibit star to star variations, 
either due to inhomogeneities at the time of for- 
mation or variations in the evolutionary impact of 
post-main-sequence internal adjustments. If the 
source of the anomalous indices in NGC 6253 is 
an enhanced distribution of a-clcments, it is com- 
mon to all the stars in the cluster and must be 
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the product of galactic chemical evolution. To our 
knowledge, NGC 6253 is the first open cluster to 
exhibit potential evidence for a non-scaled-solar 
element distribution. 

How does this solution mesh with observations 
of the Galactic disk interior to the sun? Two ob- 
servational piccc!s of evidence lend some support 
to this option, though the significance of each re- 
mains indeterminate either because the result is 
too recent to be considered established or because 
the conclusions remain controversial. One piece 
of corroborative evidence is supplied by spectro- 
scopic analysis of Ccphcids interior to the solar 
circle by Andrievsky et al. (2002a), indicating 
that Cepheids that lie more than 1.5 kpc closer 
to the galactic center than the sun have [Fc/H] 
between -1-0.2 and -1-0.4, significantly higher than 
those found in the solar neighborhood. Prom sim- 
ple chemical evolution models, the dramatically 
higher [Fe/H] should be typical of recently formed 
stars within this galactocentric radius. It must be 
emphasized that though the error bars on [Fe/H] 
are small, so is the sample. If this inner disconti- 
nuity is real, it needs to be confirmed with a larger 
sample extending closer to the galactic center and 
an explanation for its apparent absence from HII 
region analyses (Henry & Worthey 1999) would 
be useful. 

Another piece of the puzzle is supplied by an- 
other well-studied, older open cluster located at a 
galactocentric distance comparable to NGC 6253 - 
NGC 6791, with a controversial history of redden- 
ing, mctallicity, and age estimation. The most re- 
cent comprehensive analysis of the cluster is given 
by Chaboyer ct al. (1999), where CMD compar- 
isons to theoretical isochrones restrict the cluster 
metaUicity to [Fe/H] = -1-0.4 ± 0.1, in excellent 
agreement with the sole high dispersion spectro- 
scopic abundance analysis of a horizontal branch 
star by Peterson & Green (1998). (For an alter- 
native view of this discussion, see Taylor (2001).) 
If the single spectroscopic result can be trusted 
as typical of the cluster, it is intriguing to note 
that the star shows enhancements relative to solar 
ratios of Mg, Si, Na, and N, but not Ca. Even 
accepting the high and enhanced metallicity for 
NGC 6791, its unique properties within the clus- 
ter sample, combined with its highly eccentric or- 
bit (Scott et al. 1995), have led some to ques- 
tion where it fits within the population mix of the 



Galaxy, thin disk, thick disk, or bulge. If the last 
alternative is correct, the cluster is not relevant 
to understanding the evolution of the disk. If our 
conjecture regarding the metallicity distribution of 
NGC 6253 is correct, a trend of high metallicity 
beginning at least 8 Gyr ago and extending to the 
present appears established for the disk interior to 
~7 kpc from the galactic center, keeping in mind 
that the current location for the clusters may not 
be representative of where they formed. 

The progress in this project would not have 
been possible without the time made available by 
the TAG and the invariably excellent support pro- 
vided by the staff at CTIO. NDL gratefully ac- 
knowledges the summer support provided by the 
Tombaugh Fund at the University of Kansas. The 
cluster basic parameters project has been helped 
by support supplied through the General Research 
Fund of the University of Kansas and from the De- 
partment of Physics and Astronomy. 
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A. Broad-Band Photometric Merger 

Broad-band CCD data are available from the three photometric surveys with UBVRI included in BR 
and SA, and BVI found in PI. SA supply a summary of the differences in the indices among the various 
studies, finding significant offsets for a number of the colors. The risks in the standard approach of averaging 
the residuals for all the stars in such comparisons are: 

(a) The residuals may be color-dependent and will not show significant variation with V if the color range 
is large at all V, though the dispersion will be increased; 

(b) The number of stars and the errors invariably increase with increasing V. Real trends among the 
residuals may remain hidden if they exist at a level comparable to the errors among the faintest stars. 

The first indications of a potential problem in the broad-band data arose during our attempts to compare 
our V data, tied to the y filter, with the broad-band sources. Using all the stars in common, SA noted 
systematic differences in V in the sense: (SA - BR) = 0.13 ± 0.10 for 437 stars between V = 11 and 21 and 
(SA - PI) = -0.07 ± 0.05 for 115 stars between V = 11 and 18. The residuals appear reasonably consistent 
with the errors quoted by the investigations since the samples tend to be dominated by the fainter range of 
the V distributions. Similar comparisons were made between our V data and the above investigations and 
it became apparent that in the case of BR, while the internal errors were the best of the three studies, a 
color dependence existed among the residuals in V. 

To avoid potential confusion caused by the inclusion of fainter stars with larger errors, the color dependence 
for each survey was tested using only stars brighter than V = 16. Table 2 summarizes the statistics for the 
comparisons, including the number of stars in each sample, the number of stars excluded for deviating by 
more than 0.1 mag from the mean of the residuals, the color relation derived in the sense (Table 1 - study) 
= [a X {B — V)] + b, the correlation coefficient, R, the dispersion in the residuals for V brighter than 16 
after correcting for offsets and color terms, and the residuals after correction for stars in common at all 
magnitudes. The residuals for the transformed F as a function of V for all three studies are shown in Fig. 
12. 

The internal errors for BR are excellent, with an increase of only about one-third in the dispersion of the 
residuals in going from the brighter stars {a = 0.016) to the entire sample {a = 0.021). In contrast, both 
SA and PI exhibit a 50 % growth between the bright stars and the entire sample and a larger dispersion in 
both cases than found for BR. It should be noted that the offset and dispersion in V found by SA in their 
comparison with PI are in excellent agreement with the relative values derived from the data in Table 2. 
Unfortunately, no two photometric studies to date produce V consistent with each other at the 0.01 mag 
level, though our system lies near the mean of the broad-band V data. 

Though no simple comparison is available between the remaining uvby colors and the broad-band fil- 
ters, the broad-band indices can be of use in refining the cluster parameters, particularly since theoretical 
isochrones are generally available on such systems. To optimize the effectiveness of such analyses, a decision 
was made to merge all the data from the published surveys. Because the photometry in all cases has been 
calibrated by transforming each magnitude from an instrumental system to the standard system using an 
offset and a color term prior to calculating the color indices, the published color indices were returned to 
individual magnitudes and the differences between the magnitudes in each filter calculated. The residuals 
were then checked for an offset and any evidence of a color term. Because it had the largest overlap with the 
two remaining studies, the photometry of BR was used as the reference source, though it was not adopted 
as the standard system. The reason is that the color term found for V carries over to an even greater 
degree in B, as summarized in Table 2. The size of the color term for B is large, but identical from the 
comparisons with SA and PI. Moreover, the size of the term is comparable to the color term derived by BR 
in transforming the instrumental B magnitudes to the standard system, a color term that is dramatically 
larger for B than any other filter, including U. 

For the standard colors, the photometry of SA has been chosen by applying the same zero-point shift 
as used in transforming V to our system to all the broad-band filters. Color terms were only included in 
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transforming the V, B, U, and R data of BR to the standard system. In all other cases, the color term was 
set to zero and a simple offset rederived from the data used to generate Table 2. The magnitudes for each 
filter were then merged, giving the photometry of BR twice the weight of the photometry of SA and PI 
individually. 
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Fig. 1. — Standard errors of the mean for the V, 
b—y, mi, ci, hk and H/3 indices as a function of V. 
Major tick marks on the y-axis are 0.2 mag apart. 
The vertical scale has been offset by 0.2 mag for 
each index for visual clarity. 

Fig. 2. — Color-magnitude diagram for stars with 
at least 2 observations each in b and y. Filled 
circles are stars with an indication of variability. 

Fig. 3.— CMD for stars within 100" of the cluster 
center. Crosses are stars with errors in b—y greater 
than 0.010 mag. 

Fig. 4. — CMD of stars in cluster core using u — y 
as the color index. Triangles are stars tagged as 
nonmcmbers blueward of the cluster mean rela- 
tion while crosses identify binaries and nonmem- 
bers redward of the mean relation. 

Fig. 5. — Traditional CMD for the stars in Fig. 4. 
Symbols have the same meaning. 

Fig. 6. — H/3 versus 6 — y for highly probable clus- 
ter members at and blueward of the turnoff; the 
three filed circles denote the positions of stars be- 
lieved to be foreground stars. Solid curve is the 
standard relation reddened by E{b — y) = 0.190. 

Fig. 7. — Relation between B — V from the merged 
photometry and 6 — t/ for stars brighter than V = 
16.5. 

Fig. 8. — Comparison of the cluster CMD to 
the overshoot GEN isochrones for (Y,Z) = (0.34, 
0.04). The reddening and apparent distance mod- 
ulus are E{B - V) = 0.26 and (m - M) = 11.6. 
Filled circles are probable core members, crosses 
arc probable binaries, triangles arc the remain- 
ing stars within the core and squares are all stars 
brighter than V = 13. 

Fig. 9. — Same as Fig. 8 for the scaled-solar 
isochrones of PAD with (Y,Z) = (0.32, 0.04). 

Fig. 10. — Same as Fig. 8 for the enhanced-a 
isochrones of PAD with (Y,Z) = (0.32, 0.04). 

Fig. 11. — Same as Fig. 8 for the enhanced-a 
isochrones of PAD with (Y,Z) = (0.39, 0.07). la- 
belfigl2 

Fig. 12. — Residuals in V between Table 1 and the 



cited reference after adjusting the published data 
for an offset and a color term, if necessary. 
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Table 2. Summary of Photometric Comparisons 



Survey Ref. a a a b at, R gr N Nj-ej 



AV = a{B -V) + h; (Table 1 -Ref) 



BR {V < 16) 


0.059 


0.004 


-0.106 0.004 0.634 0.016 


267 


7 


All 






0.001 0.021 


830 


46 


PI {V < 16) 


0.000 


0.006 


-0.021 0.006 0.006 0.019 


118 


4 


ALL 






-0.004 0.028 


228 


7 


SA {V < 16) 


0.008 


0.008 


0.038 0.009 0.104 0.021 


80 


8 


All 






0.002 0.033 


173 


24 




AB = 


a{B - 


y) -1-6, (BR - Ref); B < 16.5 






PI 


-0.179 


0.015 


0.184 0.014 -0.754 0.034 


118 


8 


SA 


-0.175 


0.009 


0.243 0.009 -0.908 0.021 


83 


10 




AU = 


a{B - 


V)+b, (BR - Ref); U < 17.0 






SA 


-0.070 


0.021 


0.022 0.019 -0.364 0.041 


76 


5 




AR = 


a{V- 


R) + b, (BR -Ref); R< 16.0 






SA 


-0.104 


0.020 


0.168 0.012 -0.464 0.031 


103 


13 




A7 = 


aiV - 


7) -1-6, (BR - Ref); I < 15.5 






PI 


0.033 


0.012 


-0.060 0.014 0.214 0.032 


155 


12 


SA 


0.022 


0.013 


-0.009 0.015 0.163 0.032 


109 


12 



